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In this work, we consider the solvability ol the Fokker-Planck equation with both time-dependent 
drift and diffusion coefficients by means of the similarity method. By the introduction of the simi- 
larity variable, the Fokker-Planck equation is reduced to an ordinary differential equation. Adopting 
the natural requirement that the probability current density vanishes at the boundary, the resulting 
ordinary differential equation turns out to be integrable, and the probability density function can 
be given in closed form. New examples of exactly solvable Fokker-Planck equations are presented, 
and their properties analyzed. 
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I. INTRODUCTION 
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The Fokker-Planck equation (FPE) is one of the basic tools which is widely used for studying the effect of fluctuations 
in macroscopic systems [1]. It has been employed in many areas: physics, chemistry, hydrology, biology, even finance, 
and others. Because of its broad applicability, it is therefore of great interest to obtain solutions of the FPE for various 
physical situations. Many methods, including analytical, approximate and numerical ones, have been developed to 
solve the FPE. 



i-^ ■ Generally, it is not easy to find analytic solutions of the FPE. Exact analytical solutions of the FPE are known 
for only a few cases, such as linear drift and constant diffusion coefficients In most cases, one can only solve 
the equation approximately, or numerically. The well-known analytic methods for solving the FPE include a change 
of variables, perturbation expansion, eigenfunction expansion, variational approach, Green's function, path integral, 
moment method, and the continued- fraction method (for a review of these methods, see eg. Symmetry methods 
have also been quite useful in solving the FPE [2[ . Two interesting analytic approximation approaches are the WKB 
analysis [l| and the normal mode analysis Q . The finite-difference [|| and finite-element methods [5[ are among some 
' useful numerical methods. Most of these methods, however, are concerned only with FPEs with time-independent 
, diffusion and drift coefficients. 
£f) ■ Solving the FPEs with time-dependent drift and/or diffusion coefficient is in general an even more difficult task. 
It is therefore not surprising that the number of papers on such kind of FPE is far less than that on the FPE with 
time-independent coefficients. Some recent works on the FPE with time-dependent diffusion coefficients appear in 
[^-Q, and works involving time-dependent drift coefficients can be found in [9rfllj]. Refs. [l2l - fl4j consider FPEs 
with both time-dependent diffusion and drift coefficients. The symmetry properties of the one-dimensional FPE with 
arbitrary coefficients of drift and diffusion are investigated in [l4 |. Such properties may in some cases allow one to 
transform the FPE into one with constant coefficients. It is proved that symmetry group of these equations can be of 
■ one, two, four or six parameters and the criteria are also obtained. 

In our previous work [Tll j . we have considered, within the framework of a perturbative approach, the similarity 
solutions of a class of FPEs which have constant diffusion coefficients and small time-dependent drift coefficients. 
Motivated by our previous work, in this paper we would like to study the solvability of the FPE with both time- 
dependent drift and diffusion coefficients by means of the similarity method. 

One advantage of the similarity method is that it allows one to reduce the order of a partial differential equation 
15]. Thus the FPE can be transformed into an ordinary differential equation which may be easier to solve. However, 
for the FPE to admit similarity solutions, it must possess proper scaling property under certain scaling transformation 
of the basic variables. This all boils down to the scaling behaviors of the drift and diffusion coefficients. To study 
these scaling behaviors is the main aim of the present work. 

The plan of the paper is as follows. Sect. II discusses the scaling form of the FPE. Then in Sect. Ill the similarity 
method is applied to reduce the FPE into an ordinary differential equation, and the exact solution of the probability 
density expressed in closed form. Some new examples of exactly solvable FPEs with similarity solutions are presented 
in Sect. IV and V. Sect. VI concludes the paper. 
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II. SCALING OF FOKKER-PLANCK EQUATION 

The general form of the FPE in (1 + l)-dimcnsion is 

dW(x,t) 



dt 



^D^(x,t) + ^D^(x,t)]w(x,t) , 



(1) 



where W(x,t) is the probability distribution function, D^\x,t) is the drift coefficient and D^ 2 \x,t) the diffusion 
coefficient. The drift coefficient represents the external force acting on the particle, while the diffusion coefficient 
accounts for the effect of fluctuation. W(x,t) as a probability distribution function should be normalized, i.e., 
/domain W( x > t) dx = 1 for t > 0. The domains we shall consider in this paper are the real line x € (-co, oo), and the 
half lines x £ [0, oo) and (— oo, 0]. 

As mentioned in the Introduction, it is difficult, if not impossible, to find exact solutions of the general FPE with 
time-dependent drift and diffusion coefficients. Here we shall be content with a more modest aim by seeking a special 
class of solutions, namely, the similarity solutions of the FPE. Such solutions are possible, provided that the FPE 
possesses certain scaling symmetry. Below we shall show that, if both the drift and diffusion coefficients assume proper 
scaling forms, then the FPE can be solved with the similarity method. 

Consider the scale transformation 



e x 



t = £°t, 



(2) 



where s, a and b are real parameters. Suppose under this transformation, the probability density function and the 
two coefficients scale as 



W(x,t) = e c W(x,t), D^(x,t) = e d D^(x,t) , D^[x,t) = e e D {2 \x,t). 
Here c, d and e are also some real parameters. Written in the transformed variables, eq.([T]) becomes 



, dW 
dt 



l - c - d ^D^{x, t) + e 2a - c - e -^D& [x, t)] W(x, t) 



(3) 



(4) 



One sees that if the scaling indices satisfy b = a — d = 2a — e, then eq.(Q| has the same functional form as eq.([T|). In 
this case, the FPE admits similarity solutions. We shall present such solutions below. 



III. SIMILARITY METHOD 



The similarity method is a very useful method for solving a partial differential equation which possesses proper 
scaling behavior. One advantage of the similarity method is to reduce the order of a partial differential equation 
through some new independent variables (called similarity variables) , which are certain combinations of the old inde- 
pendent variables such that they are scaling invariant, i.e., no appearance of parameter e, as a scaling transformation 
is performed. 

In our case, the second order FPE can be transformed into an ordinary differential equation which may be easier 
to solve. Here there is only one similarity variable z, which can be defined as 

x a 

z = — , where a = — and a . b ^ . (5) 

t a b r w 

For o,6^0, one has a ^ , oo. In what follows, we derive the scaling forms and closed form solutions of the 
probability and current density functions. 



A. Probability distribution function 

The general scaling form of the probability density function W(x,t) is W(x,t) — (t 1 / x Xl )y(z) , where Ai and #i are 
two real parameters and y(z) is a scale- invariant function under the same transformation @. From the assumed scaling 
behavior of W(x, t) in we have —c/a = X% — (8%/a). Without loss of generality and for clarity of presentation, we 
hereby assume the parameters (Ax,<5i) = (0, ac/a). This gives 



W{x,t) =t a Zy(z) 



(6) 
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The normalization of the distribution function is 

W(x,t) dx 

domain 



y{z) 



dz = 1 



(7) 



/ domain 

For the above relation to hold at all t > 0, the power of t should vanish, and so one must have c = —a, and thus 

W( Xl t)=t- a y(z). (8) 

Next we determine the scaling forms of the drift and diffusion coefficients D^ x \x,t) and D^ 2 \x,t). Following the 
same way of determining the transformation between W{x 1 t) and y(z), we can have D^(x,t) = (t s * /x X2 )pi(z). The 
scaling behavior for D^(x,t) leads to d = bS 2 — o,\ 2 . For simplicity, we set A2 = 0. This gives 8 2 = a — 1, where 
the relations between the scaling exponents, namely, b = a — d and a = a/b, have been used. A similar procedure is 
applied to determine D^ 2 \x,t) in terms of p 2 (z). To summarize, the scaling forms of the coefficients are 



D {1 \x,t) = t a - l Pl {z) , D^{x,t) = t 2 
With eqs. ©, © and ©, the FPE is reduced to 



1 P2(Z) 



p 2 (z)y"{z) + 2p' 2 (z) - Pl (z) + az y'{z)+ p 2 \z) - p\(z) + a y(z) = 



(9) 



(10) 



where the prime denotes the derivative with respect to z. Thus it is seen that the solvability of the FPE, eq. ([1}, 
under the similarity method depends solely on that of eq. (fTOj) . An easy way to find exact solutions of eq. (|10p is to 
relate it to either the hypergeometric equation or the confluent hypergeometric equation. This requires that p\(z) be 
a linear function of z, and p 2 (z) a (linear) quadratic function of z for the (confluent) hypergeometric case. 

However, it turns out that pi(z) and p 2 [z) need not be so restricted in order to make eq. (|10|) exactly solvable. 
This is because eq. (fT0| is exactly integrable. Integrating it once, we get 



p 2 {z)y'{z) + \p! 2 {z) - pi(z) + az] y(z) = C, 



(11) 



where C is an integration constant. 

To fix the constant C, we consider the boundary conditions of the probability density W(x, t) and the associated 
probability current density J(x,t). 



B. Relation between probability density and probability current density 

From the continuity equation 

we have 

Or 

J(x, t) = D {1) (x, t) W(x, t) - — L> (2) (x, t) W(x, t) 
Using eqs. ©, © and (0, we get 

J(x,t) = t- 1 [( Pl (z) - p' 2 {z)) y(z) - P2 {z)y\z)] . 
From eq. dill) , we can reduce the above equation to 

J(x,t) = - [azy(z) - C] = - [axW{x,t) - C] , 



(12) 



(13) 



(14) 



(15) 



where W(x,t) = t~ a y(z) and z — x/t a have been used in obtaining the second expression. 

For the domains which we are interested in, i.e., the whole line and the half lines, the natural boundary conditions 
are 



J{x, t) (boundary — 0; Zy(z^ | boundary — xW (x , £) | boundary — 0- 



(16) 
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These conditions imply that C = 0, and that J(x, t) is proportional to W(x, t) and x, 

I (XX 

J(x,t) = -azy{z) = — W(x,t). (17) 

Eq. (ITT]) can also be obtained by scaling consideration as follows. Under the scale transformation, the scaling 
behavior of J(x,t) takes the form J(x,t) = e h J(x,t), where h = —b is determined from expression (|13[) with the 
help of ([2]), (|3]), the relation b = a — d = 2a — e and c = —a. The general scaling form of the probability current is 
J(x,t) — (t Ss / x Xs ) T,(z) . Without loss of generality, we choose the set of parameters (A 3 ,<5 3 ) = (0,-1). This gives 

J(x,t) =t- 1 E(z) . (18) 

Inserting the similarity solution W(x,t) — t~ a y(z) and (fTSf into eq. (p~2|) . one finds ay(z) + azy'(z) — S'(z). This 
can be integrated to give E( z ) = J [ay(z) + azy'(z)] dz + constant = azy(z) + constant. Adopting the conditions (|16|) 
then gives eq. (fT7| . 

C. Analytic expression of similarity solution W(x,t) 

With C = 0, eq. (JTTJ) is reduced to 

p 2 {z)y\z) + [p' 2 (z) - Pl (z) + az] y(z) = 0. (19) 

This equation is easily solved to give 

Pi( z ) - P2O) - az " 



( f z 

y(z) oc exp I / dz 1 



P2{z) 

Thus the probability density function W(x,t) is given by 



P2 (z)^0. (20) 



W(x,t) = At- a exp\^J dzf(z) 

n \ - - P2( z ) - az / s / n / 01 n 

f{z) = r~\ ' p^y-i ( 21 ) 

Pi\z) 

where A is the normalization constant. It is interesting to see that the similarity solution of the FPE can be given in 
such an analytic closed form. Exact similarity solutions of the FPE can be obtained as long as p\{z) and p-iiz) are such 
that the function f(z) in eq. (f2~Tj) is an integrable function and the resulted W(x,t) is normalizable. Equivalently, 
for any integrable function f(z) such that W(x,t) is normalizable, if one can find a function p 2 {z) {p\{z) is then 
determined by f(z) and p2(z)), then one obtains an exactly solvable FPE with similarity solution given by (12T1) . 

Our discussions so far indicate a general way to construct exactly solvable FPE with similarity solutions. We shall 
not attempt to exhaust all possibilities here. Instead, we will only present a few interesting cases. We will mainly be 
concerned with cases related to f(z) having the forms (i) f(z) = Az + B and (ii) f(z) = A + B/z, where A and B 
are real constants. The corresponding W(x,t) are 

W(x,t) oc t~ a e^ z2+Bz , 

W{x,t) oc t- a z B e Az , (22) 

respectively. Possible choices of p\{z) and piiz) that give f(z) in form (i) and (ii) are, respectively, (i) pi(z) = piz+p.2, 
P2{z) = P4, and (ii) p\{z) — p\z + P2(z) — pz z - We will discuss these two cases in Sect. IV and V. 



IV. FPE WITH pi(z) = mz + p 2 AND p 2 (z) = fi 4 



Let us take p\(z) and P2(z) as 

Pl(z) = piZ + p, 2 , P2(z) = P4 , 



(23) 
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where /ii, /12 and [14 are real constants. This choice of pi(z) and /?2(z) generate the following drift and diffusion 
coefficients: 

jut 30 - 1 . (24) 



(25) 

/ii = a 

We shall discuss these two cases separately. 

A. fii ^ a 

For this case, the normalized solution, from eq. (|2ip . are 

W(*,t) = xl^h^A ~ ^(x-J^-Y) , (26) 
v ' V 27r/i 4 t 2Q I 2 f i i t 2a \ a -Hi J J 

where either (/14 > 0, /ii < a) or (/14 < 0, /ii > a) must be satisfied. The well-known diffusion equation is in this 
class with a = 1/2, ji\ = ji^ = and //4 > 0. 

The solution (|2l)|) is a Gaussian (or normal) distribution in x. The full width at half maximum (FWHM) of the 
solution is related to the parameters a, fit, [i\ and time t as C\ \J {^t 2a ) / '(a — ^1), where C\ — 2\/21n2. Hence, 
from eq. (|24)) . one can see that the FWHM is affected by the coefficient /ii of the term x/t of the drift coefficient 
and the coefficient /i4 of the term i 2a_1 of the diffusion coefficient respectively. Also, it can be seen that the location 
of the peak of the solution, x — (/i2i Q )/(a — /ii), depends on the parameters a, fj,i, [12 and time t. From eq. (|24p . 
it is seen that the location of the peak of the probability density can only be influenced by the drift coefficient. 
The parameter ^2 in the drift coefficient plays an important role in the determination of the location of the peak. 
For instance, if /12 = 0, then the peak stays at the origin and will not change with time. The value of the peak is 
yja — Mi)/ (27r^4t 2Q ), and is dependent on the parameters a, fi\, [14 an d time t. Thus it is affected by both the drift 
and diffusion coefficients. 

One notes that at fixed time t, W(x, t) with (12 and — 11% are the mirror images to each other with respect to the 
y-axis as the rest of the parameters are kept fixed. Furthermore, when the parameters (2a — — (14) take the place 
of (fj,i,/j,4) with the rest unchanged, the corresponding two W(x,t)'s are the mirror images to each other. One then 
finds that W(x,i) is invariant as the parameters (/ii, ^2, /X4) are replaced by (2a — fii,—(i2, — Hi)- 

As shown previously, the FWHM and the peak value are both related to (/nt 201 ) / (a — /ii). The parameter a can 
be either positive or negative. When a > 0, the FWHM of the solution (|26|) is getting broader, the peak is turning 
smaller and is moving away from the origin with time. On the other hand, if a < 0, then the FWHM of the solution 
will shrink, the peak value will become higher and move toward the origin as time elapses. The situation where the 
FWHM is neither expanding nor contracting is impossible because a/0. 

FIG. [1] shows the evolutions of the probability density in solution (|2T)1) for a set of the parameters, where /ii 7^ a. 

For the set of parameters taken in FIG. [TJ the drift and diffusion coefficient are D^'(x,t) — x/(2t) + 1 and 
D^(x,t) — t, respectively. From the figures above and the analysis in Sect. IV. A, one sees that the peak of the 
probability distribution W(x,t) moves to the right as time increases. This is due to the presence of the drift force. 
Furthermore, owing to the presence of /ii and [14 in the drift and diffusion coefficient, the FWHM of W(x, t) is getting 
wider, and the peak value is getting smaller. 

FIG. [T] shows the evolution of the solution of the FPE with a set of time-dependent drift and diffusion coefficient. 
In subsections C and D, we shall study the difference in behavior of the solutions as one of the coefficients changes 
from being time-independent to being time-dependent, while the other coefficient is being kept fixed. 



D^(x,t)=^j + ^t a ~ 1 , D^(x,t) = 
From eq. (f2Tj) . the function y(z) is 

exp {^ t^ 1 ~ °^ ^ + ^ 2Z 



y(z) (x 



exp 



{%*}> 



B. 



/'1 



In this case, the normalized solution is 



W(x,t) 



M2 
fi 4 t c 



expi 



(27) 
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where it is valid in x > for (/^//m) < 0; x < for {^2/ ' Hi) > 0. 

Solution (j2"7| is the time-dependent solution of the FPE with the time-dependent drift coefficient D^{x,t) — 
ax/t + H2t al and diffusion coefficient D^ 2 \x,t) = ^t 2 ^ 1 . It possesses certain symmetry properties. For instance, 
when the ratio of H2I Hi is taken to be the same for different sets of (112, ^4), the solution (1271) is invariant. Also, when 
(/j,2, ha) are replaced with (— H2, Hi) or (H2,— Hi), the corresponding W(x,t)'s are the mirror images to each other 
with the other parameters fixed. 

In the following, one example is given to illustrate the solution (|27l) and the probability current J(x, t) of solu- 
tion (|2"7) . 

FIG. [2] displays the evolutions of the probability density in solution (|27| for a set of the parameters, where Hi = a - 
FIG. [3] shows the probability current J(x, t) corresponding to the solution (|2"T|) with the same set of the parameters 
as in FIG. H 

The time evolution of W{x, t) in FIG. [2] demonstrates the decrease of the concentration near the origin and the 
spread to the area away from the origin with time. 

From (|17p and (|27|) . one can find that the maximum probability current J{x,t) occurs at x = —fi4t a / fi2, which 
moves with time. 

The plots in FIG. [5] are the same for an infinite set of (h2,H4) with H2I Hi = —3 with other parameters remaining 
the same. It means that FPEs with drift coefficient D^(x, t) — 3x/t — 3Hit 2 and diffusion coefficient (x, t) = Hit 5 
have the same solutions for arbitrary Hi 7^ 0. 

C. Examples with the same drift but different diffusion coefficients 

Let us take H2 = 0, and consider two different values of a, namely, a — 1/2 and a = I. The case with a = 1/2 
corresponds to D^(x,t) = Hi^/t and D^(x,t) = /14, with the probability distribution 

^VW-H^M- < 28 » 

The other case, a = 1, leads to D^(x,t) = D^(x,t) = /i^t and 

= S^H^fM- < 29 » 

These two examples have the same drift force, and differ only in the diffusion coefficients: constant in the first case, 
and linear in time in the second. One can thus study how the time-dependent diffusion modifies the behavior of the 
system with constant diffusion. 

As H2 — 0, the peak of the probability distribution will not move with time. Comparing eqs. (|2"51) and (|2T))) , one sees 
that the FWHM is changed from C\ \J (2/^) /(I — 2^i) to C\ \J (^4£ 2 )/(I — ^1), and the value of the peak is changed 
from y(T— 2/ii)/(47r/i4i) to y/JT— /ii) / (2tt /i^t 2 ) . The two distributions coincide only at time t c = 2(1 — /ii)/(I — 2/zi). 

The evolution of the solutions (|28[) and (|29|) are plotted in FIG. 0] and FIG. [5] respectively with two sets of the 
parameters, which are the same except a. 

Because of \i2 = 0, the peaks in both of the figures do not move as expected. It can be seen that, before time 
t c = 3, W(x,t) in FIG. |4]has a larger FWHM and a smaller peak value than W(x,t) in FIG. [5j and after t c = 3, 
the situation reverses. That is because, when < t < 1, the diffusion coefficient /Z4 causes more diffusion than H^t 
does. Hence it leads the peak value of W(x, t) in FIG. |4]to be smaller than that in FIG. [5j However, when t > 1, the 
time-dependent diffusion coefficient fat starts to have more influence than the constant diffusion coefficient /14 does. 
When t = 3, both W(x,t) has the same value for all x. After t = 3, the peak value of W(x,t) in FIG. [5] is always 
smaller than W(x, t) in FIG. |U 

D. Examples with the same diffusion but different drift coefficients 

When the set of parameters a = 1/2, and fii = H2 = is taken, one has D^(x,t) = 0, D^(x,t) = /i 4 and 
probability distribution 



W{x,t) = 
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The other set of parameters a = 1/2 and p\ = leads to D^(x,t) = p%/\/i, D^ 2 \x,t) = p^ and the probability 
distribution 

WM - - ^rt( x - 2 ^) 2 } (31) 

In this case, (j30|) is the solution of the diffusion equation with a constant diffusion coefficient p^. The solution (|30|) 
demonstrates only the effect of the diffusion. Eq. (f5T|) gives the solution with additional drift p%/^/t. The probability 
distribution ([3"Tj) moves with time but keeps its FWHM at the value of 4-^/(ln2) yx^f , the same as the solution ([30]) . 
It is easy to see in ([31]) that when P2 takes the positive (negative) value, W(x, t) will move in the +x (—x) direction 
as time elapses, and the peak value a/1/(47t/X4 t) gets smaller with time. 

The evolution of these solutions ([50)) and ([3"Tj) are plotted in FIG. [5] and FIG. [7] respectively with two sets of the 
parameters, which are the same except p2- 



V. FPE WITH pi(z) = p x z + P2 AND p 2 {z) 



M32 



The next interesting exactly solvable example we shall discuss is an FPE with pi{z) — \i\z + P2 and p2(z) — p%z. 
The corresponding drift and diffusion coefficients are 



D (1 \x,t)=p 1 -+p 2 t a - 1 , D^(x,t) =p 3 xt° 



(32) 



Eq. ([2T|) is integrable and gives 



W(x,t) 





t'2 


a— All 




H3 t a 





t-) 



X ^3 



_1 exp I 



a - fJg 
t a 



(33) 



The form of W(x, t) implies that the domain of x is defined only on half-line. For definiteness we shall take x € [0, 00). 
Normalizability of W(x, t) then requires 



a — pi 
Ms 



>0, 



M3 



> 1. 



(34) 



It is seen that the distribution W(x, t) is invariant under the changes (p2, fJ-3, x) and (—^2,-^3, —a;) (other param- 
eters being the same). This means that the distributions W(x,t) with the parameters (^2,^3) an d (—^2,-^3) are 
mirror images of each other with respect to the y-axis. It is also worth noticing that for a given a, one will have the 
same solutions for different drift and diffusion coefficients as long as the quantity (pi — o)/^3 remains unchanged. 
Furthermore, the system is invariant under the transformations a — > —a, p\ — > p\ — 2a, t — > 1/t, and P2, P3 and /14 
unchanged. Hence the distribution with (a, p\) at time t is the same as the distribution with (— a, p,\ — 2a) at time 
1/t. 

Solution (|33[) with P2 — ps presents an interesting stochastic process. In this case, W(x, t) becomes the exponential 
function, whose peak is always located at the origin. Its peak value is \{p\ — a) / (p^,t a )\, which is dependent on 
the parameters a, p\ /13 and time t, and hence is affected by both the drift and diffusion coefficients. The peak 
at x = is increasing (decreasing) as t increases for a < (a > 0). That means, by an appropriate choice of the 
drift and diffusion parameters, one can have a situation where the probability function is accumulating at the origin. 
For such situation, the effect of the drift force is stronger than that of the diffusion, causing the distribution to be 
pushed toward the origin as time elapses. An example of such situation is depicted in FIG. [51 which demonstrates 
the evolution of solution (|33l) with a — — 2 and P2 — p%. 
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VI. SUMMARY AND DISCUSSIONS 

The FPE is an important equation in many different areas. Analytic solutions of the FPEs with both time-dependent 
drift and diffusion coefficients are generally difficult to obtain. In this paper, we have presented a general way to 
construct exact similarity solutions of the FPE. Such similarity solutions exist when the FPE possesses proper scaling 
behavior. 

The similarity method makes use of the scaling-invariant property of the FPE. By the introduction of the similarity 
variable, the FPE can be reduced to an ordinary differential equation, which may be easier to solve. The general 
expression of the ordinary differential equation corresponding to the FPE with time-dependent drift and diffusion 
coefficients is given in this paper. It is interesting to find, by the natural requirement that the probability current 
density vanishes at the boundary, that the resulting ordinary differential equation is integrable, and the probability 
density function can be given in closed form. We have presented several new examples of exactly solvable Fokker- 
Planck equations with time-dependent coefficients. Symmetry properties of the solutions are also discussed. 

Of course there are many other possibilities not presented here. These can be worked out easily by the reader 
following the method presented in this paper. 

Now we would like to briefly discuss whether there exist FPEs with similarity solutions which are related to the 
so-called quasi-exactly solvable (QES) equations. QES systems give rise to a new class of spectral problems for which 
it is possible to determine analytically a part of the spectrum but not the whole spectrum [l6l - |2"o| . Thus QES systems 
are intermediate to the exactly solvable systems and the non-solvable ones. One-dimensional QES systems related to 
the sl(2) Lie algebra have been classified into ten classes in [l7j . For example, the Class I QES model is given by the 
equation (l7l 120] 

[iz 2 y"(z) - [2az 2 - (2b + n)z - 2c] y'(z) + 

Note that the parameter N appears in the z-dependent term in the coefficient of y{z). For each N = 0, 1, 2, . . . , there 
are only N + 1 exactly known solutions, and hence it is only QES. 

Now one is interested to know whether there exists an FPE whose similarity solution is related to the sZ(2)-based 
QES models. To answer this question, let us note that if a given second order differential equation of the form 
P(z)y" + Q(z)y' + R(z)y — can be cast in the form in eq. (fTUf. then one must have Q'(z) — P"(z) + R(z). It is easy 
to see that eq. (|35p does not satisfy this requirement, and thus there is no FPE which can be related to the Class I 
QES model. In fact, it can be checked that this is also true for the other classes in (l7| . 
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FIG. 1: Plot of W(x,t) versus x for solution ((26| with a = 1, /xi = 1/2, /x 2 = 1, fli = 1 and time t = 1.0, 2.0, 3.0, 4.0. 
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FIG. 2: Plot of W(x,t) versus x for solution (127[) with a = 3, /ii = 3, fi2 = —6, ^14 = 2 and time t = 1.0, 1.3, 1.6, 1.9. 




FIG. 3: Plot of J(x, t) versus x for solution (|27[1 with the same set of parameters as in FIG. [2] 
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FIG. 4: Plot of W(x,t) versus x for solution ((28jl with a = 1/2, /ii = 1/4, ju 2 = 0, ^ 4 = 1, and time t = 1.0, 3.0, 5.0, 7.0. 
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FIG. 5: Plot of W(x,t) versus as for solution with a = 1, /xi = 1/4, /x 2 = 0, ^ 4 = 1, and time t = 1.0, 3.0, 5.0, 7.0 
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FIG. 6: Plot of W(x,t) versus x for solution {30J with a = 1/2, fii = /i 2 = 0, ^ 4 = 1, and time t = 1.0, 3.0, 5.0, 7.0. 
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FIG. 7: Plot of W(x,i) versus x for solution ([3T| with a = 1/2, fj,i = 0, ^ 2 = 1, m = 1, and time t = 1.0, 3.0, 5.0, 7.0. 




FIG. 8: Plot of W(x,t) versus x for solution (|33|) with a = 



-2, Mi = 



-3, ^2 = A*3 = V 2 , and time 1 = 0-5, 0.8, 1.1, 1.4. 



